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I.   INTRODUCTION 

In  recent  years  the  study  of  non-nuclear  vulnerability  and 
survivability  of  aircraft  flying  in  a  hostile  arena  has  become  increas- 
ingly important.  Increased  aircraft  sophistification  has  brought  about 
increased  complexity  of  aircraft  systems.  Better  surface-to-air 
defenses,  high  weapon  system  costs  and  decreasing  defense  budgets  have 
made  survivability  analysis  and  design  absolutely  necessary. 

Aircraft  fuel  cells  are  very  susceptible  to  damage  because  they 
have  the  largest  surface  area  and  volume  of  all  vulnerable  aircraft 
components.  Ballistic  threats  have  caused  extensive  fuel  cell  structural 
damage  which  can  result  in  fuel  starvation,  fires  or  explosion.  Major 
fuel  cell  structural  damage  cannot  be  tolerated  in  the  case  of  "wet 
wing"  aircraft  for  obvious  reasons.  It  is  for  these  reasons  that  the 
study  of  the  failure  mechanisms  of  fuel  cells  must  be  understood  so 
that  survivable  designs  may  be  evolved. 

Projectiles  which  penetrate  fluid  filled  cells  cause  damage  many 
times  more  severe  than  that  incurred  by  impact  with  an  empty  cell.  This 
interaction  of  fluid  filled  cells  with  projectiles  is  called  hydraulic 
ram.  The  exact  nature  and  extent  of  damage  can  be  related  to  projectile 
energy  loss  in  each  of  five  phases  of  hydraulic  ram.  These  phases  are 
usually  identified  by  the  dominant  energy  transfer  mechanism  presented. 

The  first  phase  involves  projectile  penetration  through  the  entry 
wall  material.  During  this  entry  phase  projectile  energy  is  lost  by 
"punching"  a  hole  in  the  entry  wall  material.  If  the  projectile  is  not 
tumbling  and  strikes  normal  to  the  entry  wall  the  damage  size  is 
approximately  equal  to  the  maximum  cross  sectional  area  of  the  threat. 
Tumbling  and  incident  strikes  can  produce  larger  damage  sizes  and 


correspondingly  higher  penetration  energy  loss. 

If  sufficient  projectile  energy  is  available  the  threat  protrudes 
from  the  hack  side  of  the  entry  wall  and  subsequent  fluid  displacement 
produces  a  strong  hemispherical  shock  wave  which  propagates  into  the 
tank  fluid.  This  shock  phase  of  hydraulic  ram  generates  intense 
pressure  loadings  on  the  entry  wall  but  because  of  spatial  attentuation 
these  loadings  are  localized  to  the  immediate  area  surrounding  the 
impact  point.  If  the  wall  is  thin,  and  the  threat  kinetic  energy  is 
large,  or  if  during  the  penetration  phase  the  entry  wall  has  been 
severly  cracked,  shock  phase  pressure  loadings  can  cause  extensive 
damage  to  the  entry  wall  area. 

As  the  projectile  travels  through  the  tank  fluid  a  cylindrical 
cavity  is  formed.  The  cavity  is  a  result  of  a  radial  velocity  imparted 
to  the  fluid  by  the  projectile  motion.  Inside  this  cavity  are  trapped 
air,  which  entered  during  penetration,  and  fluid  vapor.  Projectile 
form  drag  dissipates  projectile  energy  in  this  drag  phase  of  hydraulic 
ram. 

If  the  fuel  cell  has  a  depth  which  is  less  than  the  distance  re- 
quired for  the  drag  phase  to  dissipate  all  projectile  energy  (for 
example  the  case  of  a  projectile  penetrating  a  "wet  wing"  from  the  lower 
surface)  an  exit  phase  of  hydraulic  ram  exists.  The  exit  phase  in- 
volves exit  wall  pressure  loadings  due  to  projectile  wall  proximity  and 
due  to  fluid  motion  toward  the  wall  that  has  been  established  by  the 
earlier  phases.  These  loadings  usually  cause  extensive,  large  scale 
damage  to  the  exit  wall. 

The  last  phase  is  produced  when  the  drag  phase  cavity  growth  is 
halted  by  fluid  pressures  due  to  the  tank  boundaries.  The  cylindrical 


cavity  changes  to  a  spherical  shape  and  collapses.  Collapse  compresses 
the  trapped  cavity  vapors  to  a  high  pressure  resulting  in  a  re-expansion 
of  the  cavity.  Expansion  and  collapse  is  repeated  several  times  before 
the  stored  energy  is  dissipated.  The  resulting  fluid  pressures  can  cause 
large  side  wall  loadings  which  can  further  damage  the  cell.  These 
oscillatory  cavity  phase  pressures  can  destroy  large  depth  cells.  The 
cavity  phase  is  also  responsible  for  pumping  tank  fluid  out  damaged 
areas  of  small  depth  tanks.  This  release  of  fuel,  atomized  by  wall 
cracks  and  petaling,  creates  an  extreme  fire  hazard. 

All  phases  of  hydraulic  ram  must  be  considered  when  designing  "ram 
proof"  fuel  cells.  Fuel  cell  damage  and  the  dominant  hydraulic  ram 
failure  mechanism  are  dependent  upon  the  aircraft  mission  (which  deter- 
mines the  path,  energy  level  and  size  of  the  most  probable  threat)  and 
the  geometry  of  the  fuel  cell  (size,  location,  and  construction).  For 
example,  if  the  tank  wall  construction  is  "thin"  and  the  threat  energy 
is  "high"  the  shock  phase  of  hydraulic  ram  can  cause  catastrophic  damage 
of  the  entry  wall.  Tanks  with  "small"  depth  are  very  susceptible  to 
damage  produced  during  the  exit  phase  of  hydraulic  ram  while  the  cavity 
phase  contributes  only  minor  additional  structural  damage.  Fuselage 
tanks  may  have  dimensions  sufficient  to  eliminate  the  exit  wall  phase. 
The  primary  ram  damage  component  is  the  cavity  phase  in  this  case. 
Efficient  fuel  cell  protection  must  be  tailored  to  the  dominant  ram 
components.  Since  ram  damage  is  configuration  dependent  extensive 
vulnerability  testing  will  be  required  to  substantiate  the  defeat  of  a 
given  threat. 

This  report  presents  the  results  of  experimental  studies  conducted 
at  the  Naval  Postgraduate  School  to  understand  the  basic  phenomena 


associated  with  the  first  three  phases  of  hydraul 


ic  ram. 


II.   PENETRATION  PHASE  OF  HYDRAULIC  RAM 

The  assessment  of  ram  damage  to  any  fuel  cell  can  be  related  to 
the  energy  dissipated  in  each  of  the  five  phases  of  hydraulic  ram.  For 
this  reason  it  is  necessary  to  be  able  to  predict  the  projectile  energy 
loss  during  entry  wall  penetration.  In  general  this  energy  loss  is  a 
function  of  wall  parameters  and  projectile  initial  energy,  obliquity, 
material  properties,  tumbling  and  shape.  For  the  purposes  of  supporting 
NPS  experimental  studies  of  the  shock  and  drag  phases  of  hydraulic  ram, 
it  was  necessary  to  be  able  to  predict  the  penetration  energy  require- 
ments of  untumbled,  copper- jacketed,  small  arms  lead  projectiles  with 
various  masses  and  nose  shapes  when  impacting  normal  to  7075-T6  aluminum 
tank  entry  walls  of  various  thicknesses. 

Dunn's  theory  for  hypervelocity  ballistic  impact  problems  assumes 
that  "there  exists  in  each  uncons trained  material  body  at  constant 
uniform  temperature  a  dynamic  yield  strength,  a  ,  that  is  a  constant  of 
the  material."  This  dynamic  yield  strength  when  multiplied  by  the 
resisting  cross  sectional  area  of  the  projectile  gives  an  estimate  of 
the  instantaneous  retarding  force  during  penetration.  The  projectile 
equation  of  motion 

x  +-Ia  (x)  =  0  (II-1) 

m 

may  then  be  solved  to  yield  the  time  history  of  penetration  depth  and 
projectile  velocity  assuming  that  the  dynamic  yield  stress  of  the  pro- 
jectile is  greater  than  that  of  the  entry  wall.  The  work  done  during 
deceleration  yields  an  estimate  of  the  kinetic  energy  loss  during 
penetration.  For  the  retarding  force  postulated,  the  work  done  is 
equal  to  the  product  of  the  dynamic  yield  stress  and  the  volume  of 
entry  wall  material  swept  out  by  the  penetration.  The  theory,  therefore, 


predicts  that  the  penetration  energy  loss  is  independent  of  projectile 
shape.   The  minimum  initial  projectile  velocity  V.  which  will  penetrate 
the  entry  wall  is  called  the  ballistic  limit  and  is  given  by: 


V  = 
o 


2a  At 

Z—  (II-2) 


m 
where  A  is  the  damage  area  and  t  is  the  entry  wall  thickness. 


The  exit  velocity  V  of  the  projectile  penetrating  a  given  wall 
thickness  is  given  by: 


V  =  V 
e    o 


o 


The  corresponding  percentage  energy  loss  is  given  by: 

E.  -  E    200a  At  /-  \ 

Experimental  studies  conducted  by  Dunn  suggest  that  the  correct  value 

of  a  is  approximately  fifteen  times  the  shear  modulus  of  the  entry  wall. 

These  theoretical  predictions  were  compared  to  a  series  of  experi- 

2 
ments  in  which  projectile  shape,  impact  velocity,  and  mass  were  varied. 

The  tank  entry  wall  was  simulated  by  a  7075-T6  Aluminum  plate  of  thick- 
ness 0.05,  0.09,  or  0.l6  inches.  Three  projectile  nose  shapes  ranging 
from  hemispherical  to  a  1.6  fineness  ratio  ogive  were  tested.  Figure 
II-l  shows  the  components  of  the  ballistic  range  used  for  these 
experiments.  The  rifle  mount  shown  in  Figure  II-2  was  adjustable  in 
azimuth  and  elevation  and  was  designed  to  hold  rifles  up  to  and  including 
30  caliber.  The  impact  velocity  was  measured  by  using  three  chronographs 
screens  mounted  in  front  of  the  test  plate.  Each  screen  had  a  five  volt 
D.C.  source  across  it.  When  the  screen  was  broken  by  the  passage  of  the 


bullet  the  circuitry  shown  in  Figure  II-3  shaped  the  pulse  and  started 
or  stopped  a  Monsanto  101B  1MH  counter.  Time  interval  between  screens 
was  converted  to  projectile  average  velocity  to  determine  the  initial 
velocity.  A  similar  pair  of  chronograph  screens  were  mounted  behind 
the  test  plate  to  measure  the  exit  velocity.  Figure  II-U  shows  a  typi- 
cal chronograph  station. 

Early  in  the  test  phase  it  was  experimentally  determined  that  the 
size  and  method  of  supporting  the  test  plate  had  no  significant  effect 
on  the  exit  velocity.  It  is  expected  therefore,  that  the  exit  velocity 
for  hypervelocity  impact  would  be  insensitive  to  wall-fluid  foundation 
effects  that  would  be  present  in  any  fluid  filled  fuel  cell.  Several 
rounds  were  fired  through  a  single  l8-inch  square  plate  and  the 
resulting  exit  velocities  were  identical  even  if  the  point  of  impact 
was  as  close  as  two  times  the  damage  diameter  away  from  any  previous 
damage.  All  further  shots  used  a  three  by  four  inch  test  plate  to 
simulate  the  fuel  tank  entry  wall.  Figure  II-5  shows  the  experimental 
apparatus . 

Spark  shadowgraphs  of  the  penetration  were  also  taken.  The  delay 
circuitry  shown  in  Figure  II-3  was  used  to  trigger  the  spark  source 
just  after  the  penetration.  Figure  II-6  shows  the  optical  setup.  The 
light  from  the  spark  source  was  directed  to  a  film  plate  by  a  concave 
mirror  on  the  opposite  wall.  The  collimated  light  then  passed  behind 
the  test  plate  and  onto  the  film  plate.  These  shadowgraphs  showed 
that,  except  for  minor  distortion  of  the  nose  tip,  that  these  copper- 
jacketed  lead  projectiles  show  no  appreciable  deformation  during  the 
penetration.  Figure  II-7  shows  a  typical  shadowgraph.  The  yield 
strength  of  copper  is  approximately  0.39  X  10  psi  as  compared  to 


0.25  X  10  psi  for  Aluminum.   The  theory  assumes  that  the  yield 
strength  of  the  projectile  is  significantly  larger  than  that  of  the 
entry  wall.  Since  these  data  show  that  projectile  shape  was  un- 
changed it  can  be  concluded  that  the  assumption  only  need  be  made  that 
the  yield  strengths  of  the  projectile  is  equal  to  or  greater  than  that 
of  the  entry  wall  for  hypervelocity  impact  problems. 

Figure  II-8  shows  a  comparison  of  experimental  impact  and  exit 
velocities  with  the  theory  for  a  U5  grain  22.2  caliber  projectile 
impacting  plates  of  various  thicknesses.  Projectile  nose  shape  was 
found  to  have  no  effect  on  the  results;  therefore,  no  differentiation 
between  shapes  has  been  made  on  this  figure.  It  is  obvious  that 
equations  II-2  and  II-3  predict  very  accurately  the  exit  velocity  for 
hyperballistic  impact.  Figure  II-9  shows  the  resulting  penetration 
energy  comparison.  Good  agreement  with  the  prediction  equation  II-1+ 
is  obtained.  Additional  shots  with  bullet  masses  as  large  as  63 
grains  have  shown  similar  excellent  agreement  with  theory. 

The  independence  of  penetration  energy  loss  with  projectile  shape 
makes  it  possible  to  extend  these  results  to  many  cases  of  interest  to 
the  fuel  cell  designer.  It  could  be  expected,  for  example,  that  the 
penetration  energy  of  warhead  fragments  impacting  normal  to  any  inter- 
vening structure  and  the  fuel  cell  wall  can  be  predicted  accurately 
using  this  simple  theory. 


8 


M 

-p 

§ 
I 

o 

CD 

t>0 

« 

O 
•H 
■P 

ra 

•H 

a) 

pq 


\ 


* 


FIGURE  II-2.  ;  Systen 


10 


Vcc    (5  Volts  d.c.) 


5  K 

50    K 
Helipot 

-ihO.I  Mf 
Hh  0.01  v  f 


X 


16      15      14 


74123 


0      9 


2       3       4       5        6      7       8 


i-« 


Ml- 


! 


220  pf 
10   K 


Vcc 


2N3404 


Z  2  K 


DUPLICATE 
DELAY 
UNIT 


3    K 


OUT 


SIGNAL 
CONDITIONER 


3  Chronograph 
2  Screen 


OUT 


T 

Vc  c 
FIQTIP5    n-3*       Adjustable   Time   Delay  Unit  Circuitry 


11 


FIGURE  II-U.  Ballistic  Range  Chronograph  Station 
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III.  Shock  Riase  of  Hydraulic  Ram 

In  order  to  predict  fuel  cell  damage  due  to  the  shock  phase  of 

hydraulic  ram  it  is  necessary  to  relate  the  change  in  projectile  energy 

3 
to  wall  pressure  loading.  The  results  of  several  investigators  have 

shown  that  protrusion  of  the  projectile  into  the  tank  fluid  produces  a 
hemispherical  shock  wave  which  is  centered  at  the  impact  point.  Pro- 
jectile momentum  and  path  after  entry  have  no  observable  effect  on  the 
shape  of  this  shock  wave  as  it  propagates  through  the  tank  fluid.  In 
fact,  even  impacts  with  obliquity  show  a  hemispherical  shock  front 
centered  at  the  impact  point.  It  may  be  concluded  from  these  results 
that  the  shock  characteristics  may  be  predicted  by  assuming  that  the 
shock  is  the  result  of  a  point  energy  source  release. 

To  predict  fuel  cell  internal  pressures  created  by  such  an  event 

the  computer  code  in  the  Appendix  was  developed  based  on  the  theory  of 

5 
Yurkovitch  .  The  theory  assumes  that  the  tank  walls  are  rigid  and  that 

the  shock  radius  is  proportional  to  the  0.8  power  of  time.  The  one- 
dimensional  Rankine-Hugoniot  equations  are  coupled  to  the  fluid  state 
equation  to  determine  shock  front  conditions.  The  equation  of  state  is 
based  upon  the  Tait  equation  for  isothermal  compressible  liquids  when 
fluid  changes  are  produced  by  an  adiabatic  and  isentropic  process.  The 
assumed  state  equation  was: 

2 


P  C 
ro  o 

P=  — 


(£-)*  -i 


Po 


(III-D 


where  p  is  the  density  and  C  the  speed  of  sound  in  the  undisturbed 
fluid  and  n  is  a  constant  for  a  given  fluid  (n  =  7.0  for  water,  n  ■  10.6 
for  aviation  fuel. 

Nondimensionalization  of  the  variables  was  accomplished  by  the 


18 


following  scheme: 


r 
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RB(t) 

M     =  -| 

s         C 


Y(§,M   )    .^I^i 
v''   s  p 

^o 


S 

f(5,M  )  =  J2%iL 

S         PoRs   <*> 

s     »£<*> 

where  R  is  the  shock  radius,  u  is  the  particle  velocity,  p  is  the 
s 

pressure  and  e  is  the  specific  internal  energy  of  the  fluid.  By  using 
the  strong  shock  assumptions  that  the  shock  front  pressure  and  energy 
are  much  greater  than  the  undisturbed  conditions  the  equations  of 
motion  reduce  to: 


1  -  Tx  (l-^) 


fl  =  01 


gl  *  *fi  ft  "  |  )  (HI-3) 

f « -i-2  Kn-1) 

nM 
s 
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Solving  for  the  conditions  at  the  shock  front 


Y,  = 


1   1  -  -^  (\n-D 


nM 
s 


0  -  -^  ^l11"1)  <IIlJ0 


nM 
s 


and  assuming  values  of  Y,  the  set  III-U  can  be  solved  for  M  ,  f  and 

-1-  S    _1_ 

0  respectively.  The  shock  front  fluid  parameters  are  now  known  as  a 

function  of  the  shock  Mach  number,  M  ,  but  their  distribution  behind  the 

'  s 

front  was  to  be  determined. 

The  distribution  of  the  flow  parameters  behind  the  shock  is  based 
on  the  assumption  that  the  density  profile  is  given  by 

ql(mJ 

Y  =  \l  (HI-5) 

The  exponent  q(M  )  is  determined  by  equating  the  total  mass  enclosed 
s 

by  the  shock 

3    n1    2 

poRs  2ttJo  ^a?  (ni-6) 

to  the  original  mass  occupied  by  the  enclosed  volume 


3 
2t*s  po  (III-7) 
3 

with  the  result 


Si  *F  «s  -  i  (ra-8) 


3 
Substitution  of  (ill- 5)  into  (III-8)  yields 
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q(Ms)  =»  3  (\   -1)  (HI-9) 

the  nondimensional  particle  velocity  0  was  calculated  "by  substituting 
the  density  profile  into  the  continuity  equation 

«-s>lf  ♦*$♦«?•-*$  (III"10) 

s 

2 

where  (3  is  the  decay  coefficient  given  by  R  R  /R  .  The  resulting 

s  s  s 

particle  velocity  was 

0  =  0±l  (111-11) 

The  known  particle  velocity  and  density  distributions  were  then 
substituted  into  the  momentum  equation 


y.      _  .„     r^ 


**\-&  t^2-1'  t  -  p  »i + »s  i-> 


(111-12) 


L  *1  s  J 

to  solve  for  the  pressure  distribution.  The  decay  coefficient  was 
determined  by  assuming  that  the  shock  radius  is  given  by 


R  =  Ct0*8  (HI-13) 

s 

The  shock  Mach  number  was  determined  from  energy  considerations. 

The  energy  imparted  to  the  fluid  in  the  shock  phase,  E  ,  is  given  by 

R  ? 

E  *  f  S  (e  -  e  +  £)  pl+nr  dr  (iII-lU) 

o  ■  «o       o   2  r 

Using  the  normal  strong  shock  assumptions  and  the  previously  derived 

density  and  particle  velocity  distribution  equation  III-lU  becomes 

Eo  ■  ^oCoVK;  £  ^  (V-1)  *l  +  *l  2^  \?+Z*t  (m-15) 

nM 
s 
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By  choosing  M  the  dependent  varibles  0  and  Y.  are  known  and  equation 

S  J.        X 

111-15  can  be  integrated  to  yield  2E  /R  versus  Mach  number.   Given  the 

o'  s 

shock  phase  energy  release  E  it  is  then  possible  to  determine  the  shock 
Mach  number  as  a  function  of  shock  radius.  There  is,  however,  an  anomaly 
at  a  shock  radius  of  zero.  To  avoid  a  finite  energy  in  zero  volume 
reference  6  suggests  that  the  initial  radius  be  given  as 


R  - 
o 


E 
o 


2np  C  2 
o  o  _ 


1/3 

(111-16) 


7 
A  series  of  experiments  were  conducted  to  substantiate  the 

theoretical  propagation  speed  of  the  shock  wave  in  a  water  filled  tank. 

Figure  III-l  shows  the  ballistic  range  components  used  for  these  tests. 

Projectile  impact  kinetic  energy  was  measured  using  three  chronograph 

screens  which  triggered  counters.  The  time  interval  between  screens 

was  converted  to  average  velocity  and  then  to  impact  energy.  A  fourth 

chronograph  screen  was  mounted  ahead  of  the  test  tank  entry  wall.  This 

screen  provided  a  reference  pulse  for  timing  events  during  the  shock 

phase.  The  circuity  in  Figure  II-3  was  used  to  trigger  a  pair  of  spark 

sources  at  preset  times  after  the  reference  pulse.  In  this  manner,  a 

pair  of  spark  shadowgraphs  of  the  events  during  the  shock  phase  was 

obtained.  Figure  III-2  shows  the  optical  arrangement  of  the  shadowgraph 

apparatus . 

Figure  III-3  shows  the  test  tank  and  stand.  The  fuel  cell  was 

approximately  cubical  with  17  inch  inside  dimensions.  The  front  and 

back  wall  support  plates  were  machined  one  inch  thick  aluminum.  These 

plates  had  a  15  inch  diameter  circular  cut-out.   The  simulated  fuel 

cell  entry  and  exit  walls  were  mounted  inside  these  plates  with  l6 
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equally  spaced  bolts.  The  entry  wall  was  a  .090  inch  thick  7075-T6 
aluminum  sheet  with  a  prepunched  entry  hole  as  shown  in  Figure  III-^. 
The  top  of  the  tank  was  open  and  tape  was  used  to  prevent  water  leakage 
before  impact.  The  side  and  bottom  walls  of  the  tank  were  one  inch 
thick  Plexiglass  on  which  a  reference  gird  was  scribed. 

The  experiments  were  conducted  at  three  initial  projectile 
energy  levels.  Figures  III- 5  and  III-6  show  the  22  caliber  projectile 
shapes  and  energy  parameters.  By  varying  the  preset  time  delay  shadow- 
graphs of  the  shock  radius  versus  time  were  taken.  Figure  III-7  shows 
typical  shadowgraphs.  The  times  shown  in  this  figure  are  after  wall 
impact.  The  hemispherical  shock  centered  at  the  impact  point  is  readily 
apparent.  The  projectile  is  obscured  by  light  refraction  caused  by  the 
large  density  gradients  produced  in  the  stagnation  region  of  the  flow. 
McMillian  has  shown,  using  X-ray  shadowgrams,  that  the  illuminated  cusp 
can  be  used  to  determine  the  location  of  the  projectile.  Also  apparent 
in  the  later  time  shadowgraph  is  the  cylindrical  cavity  caused  by  flow 
separation  from  the  projectile  surface.  It  should  also  be  noted  that 
these  shadowgraphs  show  that  the  shock  strength  (as  indicated  by  the 
band  width)  is  not  uniform.  The  strength  is  a  maximum  directly  ahead 
of  the  projectile  with  a  weakening  as  the  wall  is  approached.  This 
effect  was  probably  due  to  projectile  forward  motion  and  entry  wall 
motion  in  responce  to  the  shock  pressure  loadings. 

Figure  III-8  through  111-10  show  a  comparison  of  the  experimental 
results  with  the  theory.  The  data  for  an  impact  energy  of  1600  in-lb 
shows  excellent  agreement  with  the  theory  assumtion  that  the  shock 
radius  varies  with  the  0.8  power  of  time.  The  higher  energy  level  data 
show,  however,  that  the  shock  radius  varies  as  0.9  power  of  time.  This 
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may  be  due  to  the  differences  in  projectile  shape.   If  the  projectile  is 
considered  to  be  a  piston  in  the  early  stages  of  the  shock  formation, 
the  rate  of  change  of  fluid  volume  displaced  by  the  piston  determines 
the  momentum  exchange  between  the  piston  and  the  fluid,  hence  the 
shock  strength.  The  low  energy  projectile  had  a  nose  shape  that  was 
more  blunt  than  the  higher  energy  projectiles.  It  is  reasonable  to 
conclude  therefore,  that  the  greater  rate  of  change  of  fluid  volume  is 

displaced  by  the  blunt  projectile  producing  a  stronger  shock.  The  0.8 

1+ 

power  law  has  been  experimentally  verified  by  Stepka  for  the  case  of 

spherical  projectiles.  The  present  data  shows  that  this  law  may  be  ex- 
tended to  the  case  of  blunt  projectiles  or  possibly  warhead  fragments. 
The  data  also  indicates  that  a  0.9  power  law  more  closely  predicts  the 
shock  propagation  characteristics  for  projectiles  with  nose  shapes  that 
approximate  an  ogive. 

Also  shown  in  Figure  111-10  (X  symbols)  are  data  for  the  case  of 
projectile  penetration  through  a  solid  0.0U3  inch  thick  6061-T6  aluminum 
entry  wall.  These  limited  results  show  that  the  shock  propagation  speed 
is  relatively  insensitive  to  the  entry  wall  penetration  process.  The 
spark  shadowgraphs  in  Figure  III-ll  show,  however,  that  the  solid  wall 
shock  phase  wave  system  is  considerably  more  complex  than  the  prepunched 
case.  Entry  wall  motion  in  response  to  the  penetration  phase  produces 
a  weak  wave  system  with  the  shape  of  a  truncated  cone  which  intersects 
the  main  shock  tangentially.  The  base  dimension  of  this  weaker  shock 
was  found  to  be  approximately  equal  to  the  product  of  the  speed  of 
sound  in  aluminum  and  time.  The  primary  shook  radius,  however,  was 
found  to  be  approximately  equal  to  that  of  the  prepunched  wall  case. 
The  resulting  wave  interactions  completely  obscure  the  projectile  in 
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this  case  and  presumably  alter  the  wall  pressure  loadings  during  the 
shock  phase.  Figures  111-12  through  lU  show  the  damage  resulting  from 
this  type  of  test.  These  figures  show  that  although  insufficient  energy 
was  available  to  cause  entry  wall  catas tropic  failure,  considerable 
plastic  deformation  and  buckling  occured.  Since  the  exit  wall  exhibited 
no  damage,  it  is  reasonable  to  assume  that  the  entry  wall  damage  was 
due  to  shock  phase  loadings  with  minor  contributions  due  to  loadings 
produced  by  the  latter  phases  of  hydraulic  ram. 

Implicit  in  the  Yurkovitch  solution  for  shock  phase  pressure 
loadings  is  a  knowledge  of  the  energy  input  E  .  An  analytical  expression 
for  this  energy  as  a  function  of  projectile,  wall,  and  fluid  characteris- 
tics is  not  available.  In  addition,  the  motion  of  the  thin  fuel  cell 
walls  in  response  to  these  impressed  pressures  will  reduce  the  fluid 
pressure  in  the  vicinity  of  the  wall.  An  estimate  of  these  effects  is 
one  of  the  ultimate  goals  of  the  present  NPS  hydraulic  ram  research 
effort . 

If,  in  the  design  effort  to  build  shock  phase  survivable  fuel  cells, 
the  most  conservative  approach  is  assumed  (all  of  the  projectile  energy 
is  lost  and  the  walls  are  rigid)  the  solution  of  Yurkovitch  can  yield 
entry  wall  pressure  loadings  until  such  a  time  that  the  shock  wave  re- 
flects from  any  other  tank  wall.  Figures  111-15  through  111-21  show 
such  calculated  pressure  profiles  for  the  three  energy  levels  of  the 
previous  experimental  data.  Figures  111-15,  l6,  18  and  20  show  the  pre- 
dicted pressure  signals  that  would  exist  at  a  constant  radial  position 
verses  time.  These  figures  show  that  the  shock  strength  is  quickly 
attentuated  spatially.  It  is  expected,  therefore,  that  shock  phase 
damage  would  be  limited  to  the  immediate  area  surrounding  the  impact 
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point.  Figure  111-22  shows  the  shock  radius  and  time  at  which  the  pre- 
dicted shock  wave  Mach  number  is  almost  sonic  (M  =>  1.05).  The  theory 
assumes  no  counter  pressure  (p  «  p).  For  this  reason  it  would  be 
expected  that  sonic  shock  speed  would  actually  be  reached  at  an  earlier 
time  and  smaller  radius  than  shown.  However,  this  figure  provides  a 
conservative  estimate  of  the  possible  extent  of  damaging  shock  phase 
pressure  loadings  and  of  the  time  scale  involved  in  the  shock  phase 
phenomenum.  These  theoretical  predictions  for  the  time  at  which  the 
shock  is  acoustic  were  22  micro-seconds  at  1.71  inches  for  the  1600 
in-lb  energy  level,  36  micro-seconds  at  2.31  inches  for  the  7^93  in-lb 
energy  and  hk   micro-seconds  at  ^.k   inches  for  the  12323  in-lb  energy 
level. 
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FIGUEE   III-2. 


Schematic   of   Shadowgraph  Apparatus 
(Down-range  View) 
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FIGURE  III-3.  Shadowbox,  Tank  and  Stand  Installation 
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FIGURE   III-U, 


Front    Wall    Test    Specimen, 
Pre-punched   Hole 
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Igure   III- 5.      Projectile    Shapes    (L.to  R.  .  .22 
caliber,  m=5.75xl0~3  lb,  .222  caliber,  m=6.43 
x|0"3|b,  .222  coliber,  m  =  7.85xl0~3  lb) 


Caliber 

22  Long  Rifle 

.222  Remington 

.222  Remington 

Impact  Velocity,  Vo 
(ft/sec) 

1200 

2500 

2900 

Projectile  Mass,  m 
<lbm> 

5.75xl0"3 

6.43xl0"3 

7.85  xlO"3 

Impact  Energy,  E0 
(in-lb) 

1600 

7493 

12323 

Experimental    Drag 
Coefficient,  Cq 

1.20 

0.45 

l 

0.50 

Frontal  Area  of 
Projectile  After 
Traversing  Fluid,  A 
(in2) 

0.122 

0.098 

0.165 

Equivalent   Flat 
Plate   Area,  f  (=CnA) 
fin*, 

0.146 

0.044 

0.083 

igure  IIT-6.       Projectile     Parameters 
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t  =  1  usee 


t  -  24  psec 


Figure  III-7. 


Shadowgraphs  After  Projectile  Impact  Through 
a  Pre-punched  Front  Wall 
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Figure  III-ll.   Shock  Waves  Resulting  from  0.22  Caliber 
Projectile  Impact  with  a  Solid  Entry 
Wall 
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Figure  111-12 .  Solid  Front  Wail  After  Test  (6061-T6  Aluminum, 
.043  in.  thick,  E0  =  12323  in-lb) 


Figure  111-13. 


Typical  Damage  to  Solid  Front  Wall  at  Point 
Impact  (Front  view,  6061-T6  Aluminum,  .043 
thick,  En  =  12323  in-lb) 
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Figure  111-15.   Pressure  vs.  Radius  from  Impact  - 
Complete  Profile,  EQ  =  1600  in-lb 
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Figure  ITI-16.   Pressure  vs.  Radius  from  Impact, 
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Figure  TTI-20.       Pressure   vs.    Radius    from   Impact, 
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Figure  111-22.  Radius  and  Time  for  Acoustic  Shock 
Speed  in  Water  Verses  Energy  Input 
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IV.  DRAG  PHASE  OF  HYDRAULIC  RAM 

After  protrusion  from  the  entry  wall  the  projectile  transfers  its 
remaining  energy  to  the  fluid  "by  the  production  of  the  shock  phase  pre- 
viously discussed  and  by  viscous  dissipation  as  it  traverses  the  tank 
fluid.  Figure  IV-1  presents  a  qualitative  view  of  this  projectile 
energy  decay.  The  shock  phase  does  not,  in  most  cases  of  interest, 
dissipate  all  of  the  projectile  energy.  The  projectile  subsequently 
travels  through  the  tank  until  the  remaining  energy  is  dissipated  or 
until  the  exit  wall  is  reached.  During  this  latter  drag  phase  of 
hydraulic  ram  a  cylindrical  cavity  is  formed  as  a  result  of  radial 
velocity  imparted  to  the  fluid  by  the  projectile  motion.  Inside  this 
cavity  are  trapped  air  which  entered  during  the  entry  wall  penetration 
and  fluid  vapors. 

During  the  drag  phase  projectile  spin  and  angle  of  attack  combine 
to  produce  Magnus  sideforces  which  perturb  the  projectile  from  a  straight 
path.  Furthermore,  large  stagnation  pressures  in  the  flow  cause  projec- 
tile deformation  and  in  some  cases  the  projectile  does  not  remain  intact. 
It  is  also  likely  that  the  energy  imparted  to  the  fluid  by  viscous 

dissipation  during  the  shock  phase  would  modify  the  pressure  distribu- 

9 
tion  behind  the  shock.  Fuhs  has  suggested  that  a  characteristic  time 

separating  the  drag  and  shock  phases  of  hydraulic  ram  be  identified  as 
the  point  in  time  where  projectile  stagnation  pressure  just  exceeds  the 
shock  front  pressure.  This  characteristic  time  was  approximately  30 
micro-seconds  for  the  l600  in-lb  energy  level  projectile  tested.  Inter- 
actions between  the  shock  phase  and  drag  phase  pressure  loadings  would 
be  largest  for  times  close  to  this  characteristic  time. 

These  factors  make  an  analytical  estimation  of  drag  phase  pressure 


hi 


k 

loading  exceedingly  difficult.  Stepka  has  suggested  that  during  the 
drag  phase  the  projectile  energy  decay  may  be  approximated  by  solutions 
of  a  simple  equation  of  motion 

D  "   "  m  dT  a  CD  \  pv2a  (IY"1) 

In  general  projectile  deformation  would  require  both  drag  coefficient 
and  area  variation  with  time.  By  assuming  constant  projectile  area, 
drag  coefficient,  and  mass  equation  IV-1  may  be  integrated  to  yield  the 
projectile  velocity  decay  as  a  function  of  distance  traveled  through  the 
fluid 

—  =  e  (IV-2) 

0 

where  B  =  pC_A/2m  and  V  is  the  initial  velocity  of  the  drag  phase. 

Projectile  position  as  a  function  of  time  is  given  by 

x  =  |  in   (BVQt  +  1)  (IV-3) 

From  this  crude  analysis  an  estimate  of  the  rate  of  energy  transfer  to 
the  fluid  may  be  made  from: 


f-  =■  1 ^  (IV-10 

E0   [l+BV0tf 

7 
A  series  of  experiments  were  conducted  to  determine  the  drag  phase 

energy  dissipation.  Using  the  experimental  setup  discussed  in  section 
III  spark  shadowgraphs  of  the  drag  phase  were  taken  for  three  initial 
impact  energy  levels.  Figure  IV-2  shows  typical  results.  The  cylin- 
drical cavity  formation  is  clearly  shown  in  these  shadowgraphs.  The  tip 
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of  the  illuminated  cusp  marks  the  position  of  the  projectile  nose.  From 
measurements  of  projectile  position  versus  time  the  energy  decay  of  the 
drag  phase  may  be  determined.  Figure  IV-3  through  IV- 5  show  a  comparison 
of  the  experimental  and  theoretical  positions  for  the  three  energy 
levels  tested.  It  was  hoped  that  this  data  would  yield  an  estimate  of 
the  energy  loss  during  the  shock  phase,  hence  the  unknown  initial  drag 
phase  velocity  V.  defined  in  Figure  IV-1.  Data  scatter,  however,  was 
too  great  to  make  such  an  estimate  and  V-  was  subsequently  assumed  to 
be  equal  to  the  impact  velocity.  Data  scatter  was  attributed  to  varia- 
tions, from  shot  to  shot,  in  projectile  path  and  deformation.  It  was 
concluded  that  trajectory  parameters  hence  drag  phase  energy  decay  is 
very  sensitive  to  projectile  attitude  and  deformation  during  fluid 
traverse. 

Approximately  51$  of  the  22  caliber  projectile  impact  velocity  was 
lost  within  1.5  inches  of  the  front  wall  and  within  135  micro-seconds. 
The  22.2  caliber  projectile  with  an  initial  energy  of  7^+93  in-lb  lost 
only  17%   of  its  impact  velocity  at  1.5  inches  while  those  with  an  energy 
of  12,323  in-lb  lost  25$.  The  differences  in  these  results  are  attri- 
buted to  differences  in  projectile  deformation  rates.  The  lead  22 
caliber  blunt  projectile  deforms  more  rapidly  that  the  copper  jacketed 
22.2  caliber  rounds.  It  would  be  expected  that  this  round  would  lose  a 
larger  percentage  of  the  initial  energy  at  a  constant  penetration  depth 
because  of  the  resulting  larger  equivalent  flat  plate  area  (f  ■  C_A). 
The  percentage  energy  lost  at  a  constant  radius  was  also  a  function  of 
the  initial  energy  level  for  projectiles  of  the  same  construction.  In- 
creased initial  energy  produces  larger  viscous  dissipation  rates  and 
larger  deformation  rates  because  of  higher  stagnation  pressure  loadings . 
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This  effect  is  shown  clearly  by  examining  typical  projectile  shapes 
after  fluid  traverse  as  shown  in  Figure  IV-6.  The  projectiles  studied 
are  known  to  "by  dynamically  stable  when  the  longitudinal  axis  is  perpen- 
dicular to  the  flight  path  and  hence  they  do  not  tumble.  Stagnation 
pressures  generated  by  the  drag  phase  flow  are  large  enough  to  cause 
large  plastic  deformations  of  these  lead  and  copper  jacketed  projectiles. 
Since  drag  coefficient  and  area  are  changing  during  the  drag  phase  it  is 
more  convenient  to  use  an  equivalent  flat  plate  area  to  match  experimen- 
tal data  to  the  predictions  of  equation  IV-U.  These  results  indicate 
that  the  exponential  velocity  decay  can  be  used  for  a  reasonable  estimate 
of  the  drag  phase  energy  loss  if  an  equivalent  flat  plate  area  of  the 
round  in  question  can  be  determined  experimentally.  These  data  also 
indicate  that  the  drag  phase  energy  dissipation  rate  is  very  sensitive 
to  the  material  and  aerodynamic  properties  of  the  threat. 
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Figure  IV-2. 


Cavity  Growth  Resulting  from  a 
Prepunched  Hole  Impact 
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Figure  17-6. 


Typical  Deformation  of  Projectiles  After 
Coming  to  Rest  in  the  Fluid  Filled  Tank  - 
Pre-punched  Hole  Results  (L.  to  R. :   .22 
caliber,  E^  =  1600  in-lb;  .222  caliber, 
E   =  7493  in-lb;  .222  caliber,  EQ  =  12323 
in-lb) 
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V.   Recommendations  for  Further  Study 

Hydraulic  ram  studies  will  continue  in  FY75  and  the  following  list 
outlines  the  future  work: 

1.  Information  similar  to  that  presented  in  this  report  will 
he  collected  using  the  higher  energy  levels  of  30-callber 
projectiles. 

2.  The  test  cell  will  he  instrumented  with  pressure  trans- 
ducers to  determine  experimental  pressure  profiles  during 
the  shock  phase  of  hydraulic  ram  for  cells  with  rigid  walls. 
A  comparison  of  these  results  with  predictions  will  allow 
an  assessment  of  the  Yurkovitch  model  and  could  provide  a 
method  for  determining  the  shock  phase  energy  loss  E  . 

3.  Fuel  cell  pressure  loadings  will  he  determined  for  cells 
with  thin  entry  walls.  A  comparison  with  the  results  in 
item  2  should  result  in  a  determination  of  the  fluid-wall 
motion  interaction  prohlem. 

k.     The  test  cell  entry  walls  will  he  instrumented  with 
strain  gages  to  determine  wall  stresses  during  hydraulic  ram. 
These  data  will  he  compared  to  theoretical  predictions  "based 
upon  "both  the  experimental  and  the  theoretical  pressure 
loadings.  The  theoretical  predictions  will  he  ohtained  from 
the  digital  computer  code  SATANB   and  BR-1HR. 
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VIII.     APHINDEC 

COMPUTER  PROGRAM 


C   ********************************************************* 

C  *  * 

C  *  * 

C  *  THIS  PROGRAM  CALCULATES  THE  FOLLOWING                 * 

C  *  PART  I                                            * 

C  *  1.  SHOCK  RADIUS  VERSES  TIME                    * 

C  *  2.  SHOCK  MACH  NUMBER  VERSES  TIME               * 

C  *  PART  II                                           * 

C  *  1.  PRESSURE  VERSE  TIME  AND  RADIUS              * 

C  *  FROM  IMPACT  BEHIND  THE  SHOCK                * 

C  *  FOR  A  GIVEN  FLUID  AND  IMPACT  ENERGY, EO.   THIS  PROGRAM* 

C  *  ASSUMES:  SHOCK  RADIUS  IS  PROPORTIONAL  TO  TINE  TO  THE  * 

C  *  C.8  POWER  UNTIL  ACOUSTIC,  A  POWER  LAW  DENSITY  PROFILE* 

C  *  BEHIND  THE  SHOCK,  A  STRONG  SHOCK,  ADIAbATIC  ACROSS    * 

C  *  THE  SHOCK,  THE  UNDISTURBED  FLUID  FRESSURE  IS  ESSEN-   * 

C  *  TIALLY  ZERO,  SPECIFIC  INTERNAL  ENERGY  IN  LNCISTURBEC  * 

C  *  FLUID   IS  ESSENTIALLY  ZERO.                           * 

C  *  * 

C  *  * 

C  ******************************************** ************* 

C 
C 

C      SYMBOL  DESIGNATION 

C  EO=  IMPACT  ENERGY  OF  PROJECTILE  (IN-LE) 

C  CO=  SPEED  OF  SCUND  FOR  FLUIC  (FPS) 

C  RHOO=  DENSITY  IN  UNDISTURBED  FLUID  (LE/FT3) 

C  DC=  INCREMENT  OF  NON-DIMENSIONAL  RADILS,  DELTA 

C  ZHETA 

C  DELT  =  INCREMENT  OF  TIME  (MICRC-SECi 

C  ET=  ENERGY/VOL.  =  2*E0/(R**3) 

C  Fl=  NON-DIMENSIONAL  PRESSURE  AT  THE  SI-CCK  VS. 

C  SHOCK  MACH  NO. 

C  DODM  =  DERIVATIVE  CF  NON-CI  P  ENS  IONAL  PARTICLE 

C  VELOCITY  W.R.T.  SHOCK  PACh  NO.  VS  SFCCK 

C  M-ACH  NO. 

C  VS=  SHCCK  MACH  NC. 

C  X=  PERCENT  EO  LOST  AT  ENTRY  kALL 

C 

c 

DIMENSION   ET(20) ,PSI1(20) ,F1(20)  ,VS(20)  ,  DCDP(20) 
1,RXMS(20),CRS(20), ARMS ( 20 ) , DT( 20 ) ,T( 20) , RE(20), 
1RW20)  ,  ROOD  (20)  ,  XMS  (20)  ,RS  (20)  ,  PPS  I  (20)  ,  RF  (20  ) 

CIMENSION   TT(300) ,XM(300) ,  R( 300 ) ,F2( 300)  ,  PSI2 (300 )  , 
1DCDM2(300)  ,Q(300)  ,C(100)  ,F(2   ,20  )  ,P(2   ,20) 
DIMENSION  C0NRS(30C),SH0CV(300) ,TCF0R(3  00) 
CIMENSICN  F0R(20) 
DC  500  J=l,20 
F0R(J)=C.O 
500  CCNTINUE 
N  =  lfi 

READ  (5,100)  EO,CO,RH00,X,DC,DELT 
READ  (5,101)(ET(I)  ,1=1, N  ) 
READ  (5,101)(PSI1(I)  ,I=1,N) 
READ  (5,101)(F1(I  )  ,  1  =  1, N  ) 
READ  (5,101)(VS( I)  ,  1=1, N  ) 
READ  (5,101)(D0DM( I)  ,I=1,N) 
WRITE  (6,200  )  EOtCOf RhCOfXfDCt CELT 
DC  15  1=1. N 
15  WRITE  (6,201)   ET (  I  )  ,PS  II ( I  ) , F 1( I  )  , VS( I )  ,DCCM( I ) 
C 

C      PART  I 
C 
C 


61 


C      PERCENTAGE  OF  INITIAL  ENEPGV  DEPOSITED  ANC  INITIAL 

C      ShCCK  RADIUS 

C 

EG=X*EO 

RC=(E0*32,.2*144./(PH00*C0*C0*2.*3.141o)  )**0. 33333 
C 

C      INITIAL  VALUE  OF  ENERGY/VOL. 
C 

ET0=2.*EC/(RC**3.0) 
C 

C  IIWFRT    INPUT    ARRAYS    FOR    USE    WITH    INTERPOLATION 

C  ROUTINE,    PIF2 

C 

DC    12    I=ltN 

J=N-I+1 

PE(I)=ET(J) 

PPSKI  )=PSI1(J) 

RF(IJ=F1( J) 

RV(I)=VS( J) 
12    RDCD(I)=DODM(J) 
C 

C      INTERPOLATE  FOR  INITIAL  VALUES  AT  SHOCK 
C 

VS0=PIF2(ET0,RE,N,RV) 

F10=PIF2(VS0,RV,N,RF) 

PSI10=PIF2(VS0,RV,N,RPSI) 

DCDMO=P I F2 (VSO,RV,N,RDCC) 
C 

C      COMPUTE  INITIAL  TIME  FOR  DISTANCE  RO 
C 

TC=R0/(CC*VS0*12.0) 
C 

C      ASSIGN  SUBSEQUENT  VALUES  OF  SHCCK  MACH  NC .  /5T  WHICH 
C      SHCCK  RADIUS  WILL  BE  COMPUTED 
C 

OC  10   1=1, N 

IF  (VSO-VS(I))   1,2,2 

1  GC  TO  10 

2  M=N-I+1 

xms(1)=vso 

GC  TO  11 
1C  CONTINUE 

c 

C      COMPUTE  SUBSEQUENT  VALUES  OF  SHCCK  RADIUS 
C 

11  RS(1)=R0 

DO  20   J  =  2,M 

xms(j)  =  vs.(I) 

RS(J)=(2.*E0  /(ET(I))  )**0.3333 
1  =  1  +  1 
20  CONTINUE 
K  =  M-1 
C 

C      CCMPUTE  DELTA  TIME  CORRESPONDING  TC  SHOCK  fACH  NO. 
C      ANC  SHOCK  RADIUS 
C 

CC  30   1  =  1, K 
DRS(I)=  RS(I+1)-RS(I) 
ARMS(I)=2-/(XMS( I )+XMS( 1+ 1 )  ) 
3  0  0T(I)=DRS(I)*ARMS(I)/(C0*12.0) 
T(1)=T0 
DC   40   J  =  2,M 
40  T(J)=  T(J-1)+DT( J-l) 
WRITE  (6,102) 
kRITE  (6,103)   EO 
WRITE  (6,110) 
C 

C      PRINT  SHGCK  MACH  NC.  ANC  SHCCK  RACIUS  VS  TINE 
C 

DC  50   1  =  1, M 
5C  fcRITE  (6,111)   T( I  ),XMS(I) , RS ( I ) 
100  FORMAT   (6E10.0) 
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1C1  FCRMAT   (8E10.0) 

102  1-GRMAT  (1H,//,39X,  • ENERGY • , 15X ,  •FLUID1  ) 
1C3  FCRMAT  <29X,1E20.7) 

11C  FORMAT  (1H,//,39X  • TI ME • ,17X, • MACF • ,17X , ' PACIUS '  ) 
111  FCRMAT   (29X,3E20.71 
2CC  FCRMAT  (//,  6E20.7) 
2C1  FORMAT  (15X,5E20.7) 
C 
C 

C      PART  II 
C 
C 

C      INITIALIZE  VARIABLES  AND  ASSIGN  VALUES  USEC  WITH 
C      NCN-DIMENSIONAL  PRESSURE  DISTRIBUTION  FUNCTION 
C 

8eB=T(M)*l.E06+l.G5 
L  =  BBB 
TT(l)=O.G 
XM(1)=XMS(1) 
R(ll«RO 
F2(1)=F10 
PSI2(1)=PSI10 
DCDM2(1)=DCDM0 
Q(1)=3.*(PSI10-1.0) 
e=-0.25 

I*1.0/DC  +1.05 
C(l)=O.C 
WRITE  (6,301) 

WRITE    (6,302)    TT ( 1 ) ,XM ( 1 ) , R (1 ) 
WRITE    (6,303) 
OC    60         J=2,N 
€C    C(J)=C(  J-U+DC 

c 

C      CCMFUTE  NON-DIMENSICNAL  AND  CIMENSICNAL  PRESSURE  FOR 

C      INITIAL  ShCCK  RADIUS. 

C 

DC  70  J=1,I 

F(1,J)=F2(1)+PSI2(1)*(C( J)**(4(1)*2.0)-1.C)*(F2(1)/ 
1PSI2(1)-B*(F2(1)+XM(1)*DGDM2(1)  )  ) / (Q( 1) +2  .0  ) 

Pll,J)=F(l,J)*RH00*CO*C0*XM(l)*XMl)/(32.2*144.) 

PUNCH  91,R(1),C( J),P(1,J) 
70  CCNTINUE 

K=l 

JJ  =  1 

SUM=0.0 

CALL  FORCE(F0R,C,P,R, I,JJ,K) 

CC  501  J=1,I 

5LM=SUM+F0R(J) 

WRITE (6,304)  C( J ) , P ( 1 , J  ) , F(  1, J ) , FCR( J ) 
5C1  CCNTINUE 

TCFOR( 1)=SUM 
C 

C      ITERATION  TO  SOLVE  PRESSURE,  SFCCK  RADIUS,  SHCCK  MACF 
C      NO.  VS  TIME 
C 

DC  80   K=2,L 

TT(K)=TT(K-1)+DELT 

IF  (TT(K)-TO)    21,21,22 

21  XM(K)=XMU) 
R(K)*R(1 
F2(K)=F2(1) 
PSI2(K)=PSI2(1) 
DCDM2(K)=D0DM2(1) 
Q(K)=  Q(l) 

CC  55   J=1,I 
55  P(2,J)=P(1,J) 
PUNCH  92,R(K) 
GC  TO  80 

22  X=TT(K) 
C 

C      INTERPOLATION  OF  INPUT  DATA  AND  PART  I  DATA  FOR 
C      VARIABLES  IN  PRESSURE  FUNCTION 


63 


XMK)=PIF2(X,1,M,XNS) 
R(K)=PIF2(X,T,M,RS) 
Y  =  XMK) 

F2(K)=PIF2(Y,RV,N,PF) 
FSI2(K)=PIF2(Y,RV,N,RPSI) 
DCDM2(K)=PIF2(Y.PV,N,RDCD) 
g(KJ=3.*(PSI2(K)-l.C) 
WRITE   (6,301) 

WRITE   (6,302)  TT(K) ,XN(K) ,R(K) 
WRITE   (6,303) 
DC  90   J=1,I 

F(2,J)=F2(K)+PSI2(K)*(C( J)** (Q(K)+2.0)-1.0 )*(F2(K)/ 
1PSI2(K)-B*(F2(K)+XM(K)*D0DM2(K) ))/(Q(K)+2.0) 
P(2,J)=F(2,J)*RH0C*C0*C0*XM(K)*XN(K)/(32.2*144.) 
PUNCH  91,RtK) ,C(J) ,P(2, J) 

90  CONTINUE 
JJ  =  2 
S0M=0.0 

CALL  FCRCE(FCR,C,P,R, i,JJ,K) 
CC  502  J=1,I 
SLM=SUM+FQP(J) 

WRITE  (6, 3-04)  C(J),P(2,J),F(2,J),FCR(J) 
5o2  CONTINUE 

TCFCR(K)=SLM 
SC  CONTINUE 

WFITE(6,305) 

CALL  CQNVER(R,XM,TCFOR,TT,L) 

91  FORMAT  (3E15.5) 
52  FORMAT  (1E15.5) 

2C1  FORMAT  ( 1H,// f36Xt ■ TI fE» ,17Xt • MACM , 17Xf • P4CIUS'  ) 
302  FCPMAT  (29X,3E20.7) 

3  03  FORMAT ( 1H , // ,39X, • P/RS ' , 17X , 'P • , 17X , • F« ,  19X, 'FORCE' ) 
3C4  FCPMAT(29X,4E20.7) 

3C5  FORMAT (lh,//,36X, •TIME',17X,'VELC«,17X,'RACLIS',17X,'T 
10FORCE' ) 
ENC 

FUNCTION  PIF2  ( X , X  LIST, N, FL 1ST ) 
C 
C      SECOND  ORDER  INTERPOLATION 


C 


DIMENSION  XLIST  (100),  FLIST  (100) 

ELIF  (P,Q,R,S,T)  =  ( (Q-P)*IS-T)/(R-Q)+S) 

IF  (X-XLIST(N))  2,1,1 

1  I  =  N-l 
GC  TO  5 

2  IF(X-XLISTd)  )  4,4,6 

4  1  =  1 

5  K  =  1 

GC  TO  3C 

6  K  =  2 

/  DC  a  I  -  1,N 

IF  (X-XLIST(I))  9,9,8 
fi  CONTINUE 

I=N 
9  1  =  1-1 
3C  BLIF1  =  BLIF(X,XLIST(I) ,XLIST( 1+1 ) , FLI ST ( I ) ,FL 1ST ( 1+1) 
1) 

10  IF  (K-l)  11,11,12 

11  PIF2  =  BLIF1 
RETURN 

12  IFUH-2J-N)  13,13,16 
12  IF  ((I-l)-l)  15,14,14 

14  IF(ABS(XLIST(I-l)-X)-ABS(XLIST(I-»-2)-X)  )  16,15,15 

15  L  =  1+2 
Gb  TO  17 

lfc  L  =  1-1 

17  BLIF2  =  BLIF  (X , XL  1ST ( I ) , XL  1ST ( L  )  ,  FL  1ST (  I  )  ,  FL  I  ST ( L  )  ) 
PIF2  =  dLIF  (X,XLIST( 1  +  1)  ,XLIST(L)  ,BLIF1,BLIF2) 

18  PETURN 
ENC 
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SLEROUTINE  FGRCE ( FOR,C , P ,R , I , J J , K) 

c  *  * 

c  *  * 

C  *  THIS  SUBROUTINE  COMPUTES  THE  FORCE  ON  THE  ENTRY  WALL  * 

C  *  OLE  TO  SHOCK  PHASE  FOR  A  FRONT  WALL  WITH  A            * 

C  *  PRE-PUNCKEO  HOLE, I.E.  PROJECTILE  COES  NOT  INPINGE  CN  * 

C  *  FRCNT  WALL  UPON  ENTRY.   COMPUTES  FCRCE  TC  EE  ZERO  FCR* 

C  *  PRESSURE  FIELD  OVER  THE  PRE-PUNCHED  HOLE  CP  WHERE  THE* 

C  *  PRESSURE  GOES  NEGATIVE.  THE  PROGRAM  ACCONPLISHES  THIS* 

C  *  BY  A  STAIGHT  LINE  APPPCXIMATI ON  BETWEEN  PRESSURE      * 

C  *  DATA  POINTS.                                          * 

C  *  RHRS=  NGN-DIMENSIONAL  RADIUS  OF  PPE-PUNCFEC  HCLE      * 

C  *  S  =  SLOPE  OF  STRAIGHT  LINE  BETWEEN  PRESSURE  CATA       * 

C  *  POINTS                                            * 

C  *  RX  =  THE  INTERCEPT  CF  R/RS , NON-C IN  ENS IONAL  RADIUS,     * 

C  *  WHERE  PRESSURE  IS  ZERO                            * 

C  *  FCP(N)=  FORCE(LB)  FOR  AN  ANNULUS  WHOSE  OLTSIDE  RADIUS* 

C  *  IS  C(N)  AND  INSIDE  RADILS  IS  C(N-l)           * 

C  *  * 

C   *  * 

C   ****************************************** ****** ********* 

C 
C 

DIMENSION  FOR (20) ,C(100) ,P(2,20)  ,R(300) 

RHPS=.375/R(K) 

IF(RHRS-1..)  100,25,25 
ICC  IF(JJ-l)  1,1,2 

1  CC  10  N=2,I 
S=(P(1,N)-P(1,N-1) )*4. 
IF(P(1,N-D)  3,4,4 

3  IF(P(1,N))  5,5,6 

5  FCR(N)=C.D 
GC  TC  10 

6  RX=(S*C(N-1)-P(1,N-1) )/S 
IF(RHRS.GT.C(N) )  GC  TO  43 
IF(RX.GT.RHRS)  GO  TC  60 
Ph=S*(RHRS-C(N-l) J+P(1,N-1) 
FCR(N)=(.5*(P(1,N)-PF)+PH)*3.1416*(C(N)*C(N)-RHRS 

1*RHRS)*R(K)*R(K) 
GC  TO  1C 
60  FCR(N)=P(1,N)*.5*3.1416*(C(N)*C(N)-RX*RX)*R(K)*R(K) 
GO  TO  1G 

4  IF(RHRS-C(N-1))  40,40,41 

4C  FCRIN)=(.5*(P(1,N)-P(1,N-1)  )+P ( 1 , N-i ) )*3  .  14  16* (C (N ) 
1*C(N)-C(N-1)*C(N-1) )*R(K)*R(K) 
GC  TO  10 

41  IF(RHRS-C(N) )  42,43,43 

42  PMS*(RHRS-C(N-l))+PllfN-lJ 
FCR(N)=(.5*(P(1,N)-FH)+PH)*3.1416*(C(N)*C(N)-RHRS 

1*RHRS)*R(K)*R(K) 
GC  TO  10 

43  F0R(N)=C.O 
1C  CCNTINUE 

RETURN 

2  DC  20  N=2,I 
S=<P(2,N)-P<2,N-1) )*4. 
IF(P(2,N-.D)  13,14,14 

13  IF(P(2,N))  15,15,16 

15  FCR(N)=0.0 
GC  TO  2  0 

16  RX=(S*C(N-1)-P(2,N-1) )/S 
IF(RHRS.GT.C(N))  GC  TO  143 
IF(RX.GT.RHRS)  GO  TO  160 

PH  =  S*(RhRS-C(N-l)  )+P(2,N-l) 

FCR(N)=l.5*(P(2, N)-PF)+PH)*3.1416*(C(N)*C(N)-PHRS 
1*RFRS)*R(K)*R(K) 
GC  TO  2  0 
16C  FCR(N)=P(2,N)*.5*3.14l6*(C(N)*C(N)-RX*RX)*R(K)*R(K) 
GC  TO  20 

14  IF(RHRS-C(N-1) )  140,140,141 
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140  FCR(N)=< .5*(P(2,N)~P(2,N-1) )+P(2,N-l) »*3-1416*(C<N J 
1*C(N)-C(N-1)*C(N-1) >*R(K)*R(K) 

GL    TO    20 

141  IF(RHRS-CIN))  142,143,143 

142  PH=S*(RNRS-C(N-1) )*P<2,N-1) 
FGR(N)*(.5*(P(2tN)-PHI+PH)*3.1416*(C(N)*C(N)-RHRS 

l*PhRS)*P(K)*R(K) 
GC  TO  20 

143  FCR(N)=G.O 
2C  CCNTINUE 
25  RETURN 

ENC 

SL6ROUTINE  CONVER ( R ,XM ,TOFOR ,TT , L ) 
C 

C      THIS  SUBROUTINE  CONVERTS  VELOCITY  AND  RACIL'S 
C      TC  KM/SEC  AND  CM  RESPECTIVELY 


C 


CINENSICN  R(300) , XN ( 300 ) , TOFOR ( 300 ) , CONRS ( 2 CO ) , SHOC V ( 
13CC),TT(30C) 

CC  510  11=1, L 

CCNRSUI)  =  R(1I)*2.54 

SFCCV( II)=XM(  II  J *4. 9/3. 281 

WRITE (6, 520)  TT( II  )  ,SHGCV ( I  I ) , CONRS < 1 1 ) ,  TCFCRUI) 
52C  F0PMAT(29X,4E20.7) 
510  CCNTINUE 

RETURN 

END 
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